Genetic Design of Enhanced Valley Splitting towards a Spin Qubit in Silicon 
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Electronic spins in Silicon (Si) are rising con- 
tenders for qubits the logical unit of quantum 
computation [l| owing to its outstanding spin 
coherence properties and compatibility to stan- 
dard electronics |^. A remarkable limitation 
for spin quantum computing in Si hosts is the 
orbital degeneracy of this material's conduction 
band, preventing the spin-i^ states from being an 
isolated two-level system [3, 4]. So far available 
samples of Si quantum wells cladded by Ge-Si al- 
loy barriers provide relatively small valley split- 
ting (VS), with the order of 1 meV or less (5l4lO|. 
degrading the fidelity of qubits encoded in spin 
"up" and "down" states in Si. Here, based on 
an atomically resolved pseudopotential theory, we 
demonstrate that ordered Ge-Si layered barriers 
confining a Si slab can be harnessed to enhance 
the VS in the active Si region by up to one or- 
der of magnitude compared to the random alloy 
barriers adopted so far. A biolo gically inspired 
genetic-algorithm search |llL Il2| | is employed to 
identify magic Ge/Si layer sequences of the super- 
lattice barriers that isolate the electron ground 
state in a single valley composition with VS as 
large as ^9 meV. The enhanced VS is preserved 
with the reasonable inter-layer mixing between 
different species, and is interestingly "protected" 
even if some larger mixing occurs. Implementa- 
tion of the optimized layer sequences of barriers, 
within reach of modern superlattice growth tech- 
niques [13|, overcomes in a practical systematic 
way the main current limitations related to the 
orbital degeneracy, thus providing a roadmap for 
reliable spin-only quantum computing in Si. 

The qubits for quantum information processing are en- 
coded in two- level quantum systems {|0), |1)}, and can be 
realized, for example, by two spin states { | t) , I i) } of an 
electron at the conduction band edge of a semiconduc- 
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tor [3, [isj. While Si enjoys a number of advantages 
over III-V semiconductors in this respect, including long 
spin coherence lifetime (associated with its weak spin- 
orbit coupling and small content of nonzero-nuclcar-spin 
isotopes), as well as advanced fabrication know-how, its 
major drawback is the (six-fold) orbital degeneracy of its 
lowest conduction band (Fig. [1^) located close to the X 
point in the Brillouin zone. This is no longer a two-level 
system determined solely by its spin, leading to consid- 
erable leakage and decoherence driven by the energetic 
proximity among the degenerate orbitals Whereas 
this six-fold valley degeneracy in the 0/i-symmetric bulk 
Si can be partially removed by application of tensile bi- 
axial strain [l^, thus, isolating the two lowest | + z) and 
I — z) components from the rest (Fig. [Ud), the creation of 
a sufficiently large energy splitting within this Z-valley 
subspace (hereby called valley splitting (VS), see Fig. [1];) 
has proven to be a challenge for the experimental real- 
ization of Si-based spin qubits Q . 
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FIG. 1: Description of valley splitting and system ex- 
plored. Schematic representation of (a) the (6-fold degen- 
erate) A-valley state of bulk Si; (b) splitting between XY- 
valley and Z-valley by tensile biaxial strain; and (c) splitting 
of Z-valleys by the abrupt interfaces of quantum well. Note 
that the confinement barrier height for the well is the band 
offset (between Si well and barrier) of valley (b-Valley), 
higher than the band offset of conduction band minima (b- 
CBM). The symmetry of Si quantum well of A'^ monolayers 
alternates from Y)2d ^ D2/1 for A'^ odd -f-^ even respectively, 
(d) Heterostructure geometry adopted in the present study. 
A [001] -oriented Si quantum well of thickness d and the sur- 
rounding barriers are coherently strained by epitaxial growth 
over a specified substrate. Both barrier and substrate are 
composed of Si-Ge based materials. 
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The geometry of the basic physical system explored 
(Fig. [TJi) includes a Si slab ("Well") interfaced by a ma- 
terial with higher conduction band ("barrier"). The VS 
of this system depends on a multitude of degrees of free- 
dom present in the realistic device growth. The Si well of 
thickness d cladded by barrier materials of composition 
Xb is coherently strained on a substrate with the planar 
lattice constant (determined by its composition Xg). 
We focus on the substrate and barrier composed of Ge-Si 
based materials, which provide better-quality interfaces 
than oxides. The barrier can be a Ge-Si random alloy 
of composition Xb or any corresponding atomistically or- 
dered structure. We incorporate monolithically the full 
system containing up to 10^ atoms/computational-cell, 
via an atomistic pseudopotential Hamiltonian [itI . [l8l |. 
solved in a plane- wave basis for each relaxed atomic con- 
figuration, which gives directly the energies {si} and 
wavefunctions of the valley states. 

Macroscopic degrees of freedom: We start by ex- 
ploring the continuum configuration-averaged degrees of 
freedom in this system, as common in the literature @- 
[Tol . [l9l - l2]1 | , finding that whereas they do not provide a 
clear avenue to major VS enhancement, their exploration 
hints at the importance of another length scale. We con- 
sider a fixed-thickness Si well embedded in the Ge-Si alloy 
barriers with varied composition Xb, on three substrates 
with different composition Xg. For each alloy compo- 
sition Xb of barriers, we calculated 20 randomly real- 
ized atomic configurations and the averaged VS is evalu- 
ated. The solid red line in Figs[2^-c shows the calculated 
configuration-averaged VS as a function of composition 
Xb- Generally, one observes an uneventful monotonic in- 
crease of the averaged VS as the barrier becomes richer 
in Ge (see also Supplementary Fig. Sib, which shows the 
VS for a few distinct Xb). Such continuum-like effect of 
the configuration-averaged alloy barriers can be under- 
stood by the gradual change of the barrier height, i.e., 
the band offset between the valley states of Si well and 
barrier ( "b- Valley" in Fig. [TJ; and Supplementary Fig. 
Sla) [l0|, 113. Although the averaged VS (red lines in 
Figs [2^-c) shows substantial dependence on the epitaxial 
strain (also see Supplementary Fig. Sib), the variation 
of the macroscopic barrier composition Xb and substrate 
composition X^ provides limited tuning of VS. 

Atomically resolved length scale: Important clues 
emerge as to the significance of the atomically resolved 
length scale and symmetry, as indicated in a recent work 
on the intervalley splittings of PbSe [l^. In principle, 
the splitting within the Z valleys is closely related to 
the interface- induced deviation from the Oh symmetry 
of bulk Si (or the D4/1 symmetry of biaxially strained Si). 
For a Si quantum well (Fig. [TJl), the interfacial perturba- 
tion potential AV with the T)2h/^2d symmetry provides 
a coupling channel between two Z- valley states, giving a 
VS magnitude in perturbation theory of 2|(-|-z|Ay| — 
To tune VS, we can engineer the magnitude and profile 
of the perturbation potential AV by varying the atomic- 
scale structure and symmetry for the well and barriers. 
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FIG. 2: Valley splitting of fixed-thickness Si well next 
to a barrier of disordered alloy and ordered superlat- 
tice. Calculated VS (in meV) as a function of the content 
of Ge in the barrier, for a 40 MLs Si quantum well embed- 
ded in Si-Ge disordered alloy (a,b,c) and ordered superlattice 
(d,e,f) barriers, on a 0%, 20%, and 40% Ge substrate. Solid 
red lines in (a,b,c) represent the configuration-averaged VS 
for alloy barriers, and empty circles represent specific atomic 
realizations (20 for each composition). In (d,e,f), the superlat- 
tice barrier has a 80-ML period, and structural configurations 
(green crosses) are generated by the genetic algorithm search 
(see Fig. 2)). Note that the VS for superlattice barriers (d,e,f) 
is in ~one-order larger scale than that of alloy barriers (a,b,c). 



The importance of the atomic scale is revealed, for exam- 
ple, in Fig. [2k-c, where the blue circles represent the VS 
obtained by resolving distinct random realizations of site 
occupations in alloy barriers. The VS ranging from to 
an upper bound of ~1 meV is in reasonable agreement 
with experiments [1, [H-Q . We can see that the VS of Si 
can vary significantly for different atomic configurations 
of barriers at the same composition Xb. This is consistent 
with the recent calculation showing that specific atomic 
arrangements at the interface region can result in distinct 
VS (however the assumed SiaGe luzonite structure is dif- 
ficult for experimental realization) [gj]. Also, the critical 
role of atomic resolution and symmetry is apparent by 
considering a system of short-period Si-Ge superlattices 
located directly on a substrate (i.e, no active Si layer in 
Fig. [TJi), where our calculated VS reaches values as large 
as several tens of meV, although the Si-Ge superlattice 
system is not the case of interest here (but may relate to 
different qubit proposals 25[). 

Inspired by these basic insights from the atomic length 
scale, we next explore in a systematic way whether and 
how atomic degrees of freedom in the Si well, barrier 
composition and structure, and epitaxial substrate could 
raise the VS. By varying the above degrees of freedom, we 
aim at identifying the rule of how the relevant physical 
factors govern VS, and use it to seek an optimized VS. 

Effect of Si well thickness: The thickness d of the 
Si well is the first obvious parameter to tune the pertur- 
bation potential AV, and thus manipulate VS. Fig. [3^ 
shows the dependence of VS on the thickness d in mono- 
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layers (MLs) for fixed pure Ge barrier from the pseudopo- 
tential calculations. We observe an overall decay in the 
magnitude of VS as the thickness d increases, while the 
VS for d with an odd (blue circles) and even (red squares) 
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FIG. 3: Dependence of valley splitting on Si well thick- 
ness. Calculated VS (in meV) as a function of the thickness 
d (in MLs) of a Si quantum well embedded in fixed pure Ge 
barrier on pure Si substrate. The pseudopotential (PP) and 
efi'ective mass approach (EMA) results are shown in (a) and 
(b), and the data at the odd and even number of MLs are 
shown in blue and red, respectively. For comparison with 
PP data, we show the EMA results as a continuum line with 
markings at the integer multiples of the MLs spacing. Note 
that the EMA results are smaller in magnitude than the PP 
results. 

In Fig. [Sb, we show the calculated VS within the ef- 
fective mass approximation (EMA) as a function of the 
continuum thickness dcon (solid black line) , as well as the 
data sampled at odd (blue circles) and even (red squares) 
atomic MLs. We find that while the EMA results with 
continuum dcon show a much faster oscillation, clearly 
they reproduce well the existence of "independent" os- 
cillation for discrete d of odd and even MLs. Thus, we 
attribute this atomic-scale odd-even independent oscilla- 
tion to a manifestation of the aliasing effect (introduced 
by sampling a function at a rate which is not fine enough 
to capture each oscillation), rather than to a symmetry 
change (see Supplementary Information for detailed de- 
scription). This understanding underlines that to gain 
an optimized VS of Si well, a well-controlled growth of 
monolayer precision is required to reach the thickness d 
at the peak of the oscillation. 

Atomically ordered superlattices barriers: The 
substantial effect of specific atomic realization for the 
disordered alloy barriers (open circles in Figs[2^-c) stim- 
ulates us to investigate the situation where the barriers 
are composed of ordered superlattices, i.e., a repeated se- 
quence of Si and Ge layers of arbitrarily assigned widths. 
We explore the system composed of a 40 MLs Si well 
(located at an even peak of Fig. [3^) embedded in the 



superlattice barriers with a period of 80 MLs (with the 
minimum stacking unit of bilayer to comply with current 
experimental growth conditions). This gives an astro- 
nomical number (2''°) of candidate layer-stacking config- 
urations of barriers, so a direct calculation for enumera- 
tion of all the candidates is not practical. We perform an 
inverse-band-structure search calculation jTll [l^ where 
the best fitness is defined by the maximum VS, and fa- 
vorable structures are selected within a genetic algorithm 
approach. Fig. shows the evolution of fitness (VS) 
with generation (evolution step). One clearly observes 
that the VS can be effectively tuned within a wide energy 
range, from negligibly small up to '^9 meV, by varying 
Ge/Si stacking sequence of superlattice barriers. Less 
than 100 generations already identify the best individ- 
uals, which remain superior for the following hundreds 
of generations, while new individuals still emerge with 
intermediate VS values. 
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FIG. 4: Optimization of valley splitting in inverse- 
band-structure calculations and effect of inter-layer 

mixing, (a) Evolution of fitness (VS) with generation in 
an inverse-band-structure calculation for a 40 MLs Si quan- 
tum well embedded in ordered superlattice barriers with 
80 ML period, on a 40% Ge substrate. The top axis 
shows the number of atomic configurations investigated dur- 
ing the evolution, (b) Calculated VS (in meV) for the op- 
timum configuration of superlattice barrier achieved in (a) - 
Ge4Si6Ge2Si6Ge4Si4Ge4Si4 . . ., as a function of the degree of 
inter-layer mixing, r) (see text). At = there is no mixing 
and T) = 0.5 means the maximum diffusion, i.e. complete de- 
struction of Si-rich or Ge-rich pattern layer. We explore two 
cases of mixing length: IML for each side of the interface - 
total of 2ML (green) ; and 2ML for each side - total of 4ML 
(red) , defining the maximum range at the interface where the 
mixing occurs. For each rj, 10 alloy realizations are randomly 
sampled and the averaged VS is shown as a line. The fiuc- 
tuation induced by different alloy realizations is within 0.2 
meV. 

Figs [2ji-f show the achieved VS of all the atomic con- 
figurations visited by the inverse-band-structure search, 
sorted in terms of the Ge content in the barriers on three 
varied substrates. It is demonstrated that a remark- 
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able VS enhancement by a factor of 5-10 is achievable 
with ordered superlattice barriers as compared to disor- 
dered alloy barriers (Figs H^-c) for all substrates. Par- 
ticularly, comparing with the maximum VS for the dis- 
ordered alloy barriers - '^1.0 meV on all the substrates, 
the maximum VS (accompanied by the optimum con- 
figuration) for the ordered superlattice barriers reaches: 

% Ge substrate: 5.7 meV (Ge4Si4Ge2Si6Ge4Si4Ge4Si2 . . .); 

%20 Ge substrate: 7.4 meV (Ge4Si4Ge4Si2Ge4Si6Ge4Si2 . . .); 

%40 Ge substrate: 8.7 meV (Ge4Si6Ge2Si6Ge4Si4Ge4Si4 . . .). 
We find that the multilayer superlattice barriers show 
larger VS around the central region, i.e., at 40-60% Ge 
content in the barrier, different from higher Ge content 
leading to larger VS for random alloy barriers. The same 
Ge content in the superlattice barriers can lead to both 
high and low VS extremes, again emphasizing the key 
role of atomistic scale ordering in controlling VS. 

The Si/Ge4 motif: Interestingly, all the optimum 
configurations identified start the barrier sequence by 
a Ge4 sub-layer. This same "magic" thickness for the 
first Ge sub-layer is also identified in the exhaustive enu- 
meration calculations for the superlattice barriers with a 
shorter period of 16 MLs (see Supplementary Fig. S2a- 
c). Similar results are obtained for a Si well with the 
thickness of 47 MLs (located at an odd peak of Fig. [3^, 
see Supplementary Fig. S2d-f). In order to better under- 
stand this, we explore a simple case - the fixed 40 MLs 
Si well embedded in Ge„Si„ superlattice barriers with n 
~ 1, 2, 4, 8, 16, as shown in Fig. [5^. We see that the 
barrier of Ge4Si4 superlattice indeed exhibits the largest 
VS (>7 meV), whereas all other barriers (including pure 
Ge) show typically low VS (< 2 meV). This indicates 
that the starting sub-layer thinner or thicker than Ge4 
seem to equally suppress VS. We unravel the underly- 
ing origin within the EMA context. Briefly, the VS in- 
duced by a Si/Ge (ascending offset) interface has oppo- 
site sign to the Ge/Si (descending offset) interface with 
the same wavefunction. Choosing the interface positions 
to match the maxima/minima of the VS at the ascend- 
ing/descending interfaces would maximize the total VS. 
It is impossible to match the interface positions perfectly 
to the incommensurate oscillations of well-thickness de- 
pendent VS (Fig. [3|) , but the Ge4 sub-layer is the closest 
we can get to this matching. Conversely, starting with 
a Ge2 sub-layer cladding the Si well, we find a destruc- 
tive interference, in agreement with the suppressed VS 
for Si2Ge2 superlattice barrier in Fig. [5^. This engi- 
neering is analogous to that of a distributed Bragg reflec- 
tor (see Supplementary Information for detailed descrip- 
tion) . But the fact that the oscillations are incommensu- 
rate with the lattice and the strong dependence of VS on 
atomic ordering makes it impossible to analytically pre- 
dict the optimal structure. For this reason, the genetic 
selection of candidate structures is an essential ingredient 
of this work. 

Previous studies correlate the VS with the electronic 
wave-function magnitude of the Z valley at the inter- 
face (interfacial |*|2) [1, [13, [H HI] and the wave- 
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FIG. 5: Results of Ge„Si„ barriers, and relation be- 
tween valley splitting and wave-function, (a) Calcu- 
lated VS (in meV) for a 40 MLs Si quantum well embedded 
in the GenSin (n = 1, 2, 4, 8, 16) superlattice barriers on a 
40% Ge substrate. The result for pure Ge barrier (n = oo) is 
shown for comparison. (b,c) Calculated VS (in meV) shown 
as a function of (b) planar-averaged (in the XY plane) wave- 
function squared magnitude at the interface zi (4'^ («/)), and 
(c) planar-averaged wave-function squared penetration into 
the barriers [j"^ '^'^{z)Az), for a 40 MLs Si quantum well em- 
bedded in random alloy and 16-ML period superlattice bar- 
riers, on a 40% Ge substrate. The averaged wave-function 
magnitude between two Z-valley states is used for such cal- 
culations. 



function penetration into the barrier region (penetration 
of l^-p) [13, [13. In Fig. [5b andEfc, we probed the VS 
as a function of interfacial j^'p and penetration of l^'P, 
for the 40 MLs Si well cladded by alloys barriers (blue 
crosses) and superlattice barriers (red pluses). Compared 
with the alloy barriers, the stronger confinement power 
of superlattice barriers give much narrower distribution 
of both interfacial and penetration of l^^p. The 
optimum VS values for the superlattice barriers emerge 
in the region of the narrowest distribution of these two 
quantities. This is related to the sharp well/barrier in- 
terface for superlattice barriers, which could in principle 
enhance the VS [l3,[23l. 

Effect of Ge-Si intermixing in barriers: Since it 
is still a challenge to grow perfectly pure sub-layer of Si 
or Ge in superlattices due to atomic inter-diffusion [isj . 
we examine how much VS is affected by the interfacial 
mixing between Si and Ge. In particular, the inter-layer 
mixing is modeled by mapping pure Si into Sii_,,Ge,,, 
and pure Ge into Gei_,,Si,, at the interfacial first few 
layers, determined by a mixing length. The parame- 
ter rj quantifies the degree of inter-layer mixing, with 
77 = corresponding to no mixing and rj = 0.5 mean- 
ing maximum mixing, i.e., complete destruction of Si- 
rich or Ge-rich pattern within this layer. Fig. |3)d 
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shows the calculated VS as a function of rj for the 
above optimized superlattice barrier on %40 Ge substrate 
(Ge4Si6Ge2Si6Ge4Si4Ge4Si4 . . ., see Supplementary Fig. 
S3 for more ordered superlattice barriers), when two 
cases of mixing lengths [2 MLs (green) and 4 MLs (red)] 
are explored. Note that the favorable Ge4 starting sub- 
layer is only partially damaged if the mixing length is 2 
ML (1 ML towards each side of the interface), while for 
the 4ML mixing length the Ge-pure layer is totally de- 
stroyed. The non-trivial, non-monotonic behavior indi- 
cates that the intermixing may lead to the formation of a 
more complex geometry which tunes VS by affecting the 
interference pattern discussed before. This is reflected in 
a surprisingly steeper suppression of VS in the shorter 
mixing length of 2 MLs compared to the longer mixing 
length of 4 MLs for small rj. Similarly, at very large rj, 
the structure becomes a complex layering of alloys, pure 
Si and pure Ge, which might keep suppressing (the case 
of 4ML) or invert the symmetry and enhance the VS (the 
case of 2ML). In both mixing lengths, for a reasonable 
degree of mixing (77 < 0.1), the rather high VS of >6 
meV is preserved. 

Advantage of atomically ordered barriers: We 
anticipate that the choice of ordered superlattice barri- 
ers instead of random alloy barriers might mitigate many 
problems of real samples. For instance, the intrinsic non- 
deterministic nature of alloys induces disorder ranging 
from the geometry of the interface plane to the inho- 
mogeneous strain fields [l^. The leakage of electrons 
tunneling through the superlattice barrier should also be 
reduced since the electronic density inside the barrier is 
much reduced. Our atomically-resolved pseudopotential 
calculations of the Si well strained on various substrates 
and interfaced with different barriers can be used to ex- 
plore the effects on VS of both global macroscopic quanti- 
ties (strain, alloy compositions, geometric well thickness), 
as well as atomic scale effects (ordering, inter-layer mix- 
ing, even-odd independent oscillation). We identify the 
critical Si well thicknesses as well as an emerging "magic" 
motif of Ge4 starting in the ordered superlattice barrier 
that causes strong coupling between \-\-z) and \ —z) valley 
states, leading to significantly enhanced VS as large as 
^^9 meV. The predicted structure is accessible within cur- 
rent experimental fabrication capabilities. This opens the 
way to fundamental understanding of the hitherto rather 
intangible Z-valley splitting in indirect-gap semiconduc- 
tors such as Si with the possible benefit of isolating single 
electron valley state for spin-based quantum computing. 



Methods 

The structures employed to optimize VS in this work 



(Fig. [TJi) involve an active Si well with the thickness of d 
MLs, cladded on both sides by a barrier consisting of Si- 
Ge based materials, including homogeneous random al- 
loy and layer-by-layer superlattice structures. The whole 
system is coherently strained on a substrate, via mini- 
mization of atomically-resolved strain with a generalized 
valence force field method (sec Supplementary Methods 
A) [2^. To comply with what is currently accessible in 
experimental growth we used two restrictions: (i) Since 
too high Ge content in substrate is known to cause dislo- 
cations in thick Si active layers to relieve excessive strain, 
up to 40% Ge content is considered in substrate, (ii) a bi- 
layer is used as the minimum stacking unit of each specie 
(Si/Ge) for the superlattice barrier. 

The energies and wave-functions of conduction valley 
states for candidate structures are calculated "on the fly" 
with the atomistic pseudopotential method, described in 
detail in Refs. The atomistic pseudopotential 

method (overcoming the well-known Density-Functional- 
Theory limitations on electronic structure calculations), 
accompanied with a plane-wave basis set and folded- 
spectrum diagonalization, allow us to accurately calcu- 
late energy splitting of .^-valley states (at the order of 
meV or lower) for numerous candidate structures with 
economic efficiency, as described in Supplementary Meth- 
ods B. 

Effective mass approach calculations were performed 
to accompany the interpretation of the pseudopotential 
results. The effective mass calculations follow essentially 
the model presented in Refs. [iSjIlll: adapted to describe 
quantum wells in first order approximation, as described 
briefly in Supplementary Methods C. 

For Si wells embedded in laycr-by-layer superlattice 
barriers, since the search space shows a combinatorial 
burst of degrees of freedom, we em ploy the developed 
inverse-band-structure approach [TTI. Il2l Isoj. i.e., a bio- 
logically inspired (Darwinian) genetic algorithm to guide 
the electronic structure calculations, with the aim at find- 
ing the optimum configuration that gives the maximum 
VS (Supplementary Methods D). 
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